Carrier-Density-Controlled Anisotropic Spin Susceptibility of Two-Dimensional Hole 

Systems 
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We have studied quantum-well-confined holes based on the Luttinger-model description for the 
valence band of typical semiconductor materials. Even when only the lowest quasi-two-dimensional 
(quasi-2D) subband is populated, the static spin susceptibility turns out to be very different from 
the universal isotropic Lindhard-function lineshape obtained for 2D conduction-electron systems. 
We discuss how the strongly anisotropic and peculiarly density-dependent spin-related response of 
2D holes at long wavelengths should make it possible to switch between easy-axis and easy-plane 
magnetization in dilute magnetic quantum wells. 



Introduction - The recently achieved ability to make 
bulk semiconductors magnetic [H-Q opens up excit- 
ing possibilities for controlling magnetism in unconven- 
tional ways. Nanostructuring offers appealing routes to- 
wards realizing magnetic semiconducting devices based 



or wave-function en- 
I11 particular, size 



on strain-induced anisotropies [4 
gineering in heterostructures 
quantization strongly affects the intrinsic spin-3/2 de- 
gree of freedom carried by valence-band states [7[ and 
is expected to stabilize easy-axis magnetism in a two- 
dimensional (2D) hole system with magnetization direc- 
tion perpendicular to its plane [§, ,9] . Our detailed the- 
oretical study shows that the behavior of p-typc mag- 
netic quantum wells is unexpectedly rich. Valence-band 
mixing turns out to sufficiently weaken easy-axis mag- 
netism at intermediate densities, triggering a transitions 
to an easy-plane magnet. Countcrintuitively, valence- 
band mixing is also responsible for reentrant easy-axis 
magnetism occurring at even higher density. Our results 
are robust with respect to changes in quantum-well pa- 
rameters. The density-driven magnetic transitions dis- 
cussed here enable new magneto-electronic device func- 
tionalities 13, U I based on electric- field manipulation of 
the magnetization in low-dimensional systems. 

Background & Aim - Magnetism is introduced into 
intrinsically nonmagnetic semiconducting materials via 
doping with magnetic ions such as Mn, Co, Fe or 

Gd. One way to generate an exchange interaction be- 
tween any two localized magnetic moments embedded in 
a conductor is provided by the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) mechanism [T^j, which gives rise to the 
effective two-impurity spin Hamiltonian 

T4 Q 



(i) 



Here l\ denotes the ith Cartesian component of an im- 
purity spin located at position r a , and G is the exchange 
constant for the contact interaction between the spin den- 
sity of delocalized charge carriers with the impurity spins. 

The properties of the charge carriers mediating the ef- 
fective impurity-spin exchange coupling are encoded in 



the static spin-susceptibility tensor [12j Xij( r a', r (s)- In a 
homogeneous electron system with spin-rotational invari- 
ance, it has a universal isotropic form that depends only 
on the dimensionality of the system [l]|, rendering the 
impurity-spin interaction ([T]) to be of Heiscnbcrg-modcl 
type. This case applies to the 2D electron systems re- 
alized by confining carriers from the conduction band 
in semiconductor heterostructures as long as inversion 
symmetry is not broken by the crystal lattice or due to 
structuring of the sample 14H17|. 

Here we focus on the properties of 2D hole systems 
whose charge carriers have an intrinsic spin-3/2 degree of 
freedom that is strongly coupled to their orbital motion 
even when inversion symmetry is intact Q • We adopt the 
Luttinger model [l8| in axial approximation 0, HH, [2(| , 
which provides a useful description of the upper-most va- 
lence band of typical semiconductors in situations where 




FIG. 1. Lowest three (each of them doubly degenerate) sub- 
bands of a two-dimensional hole system realized by a symmet- 
ric hard-wall quantum-well confinement of width d. Disper- 
sions are calculated based on the four-band Luttinger-model 
description of bulk valence-band states in axial approxima- 
tion, using band-structure parameters applicable to GaAs 
confined in [001] direction, and Eq = -7r 2 ft 2 7i/(2m d 2 ). 
Gray lines indicate the range of energies and wave vectors 
for which only the lowest subband is occupied. This is the 
regime we focus on in this work. 
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its couplings to the conduction band and split-off valence 
band are irrelevant. Subband k-dot-p theory [21, 22 1 is 
employed to obtain the lowest quasi-2D hole subbands for 
a symmetric hard-wall confinement characterized by its 
spatial width d. See Fig. Q] We consider the case where 
only the lowest 2D subband is occupied and calculate the 
spin susceptibility. Based on this result, various magnetic 
phases are shown to emerge in magnetically doped hole 
quantum wells depending on the carrier density. 

Spin susceptibility of a quasi-2D system - The spin 
susceptibility of 2D charge carriers is given by pjl] 

Xij{R,z;~R',z') = 

~iJo dte ~ vt tt 5 ^'*;*)'^''^ )])- ( 2 ) 

Here R and z are components of the position vector in 
the 2D (xy) plane and in the perpendicular (growth) 
direction, respectively. The spin density operator (in 
units of h) is defined in terms of field operators "J, 
and the Cartesian components Jj of the charge 
carriers' intrinsic angular momentum as Sj(R, z) = 
*t(R,z) Jj *(R,z). The field operators can be ex- 
pressed in terms of operators associated with general 
cigcnstates [labelled by band index n and in-plane wave 
vector k = (k x ,k y )] of the non-interacting Hamiltonian 



*(B, *) = £„/ 



d A k 



inv.(z) c„ k . The (normal- 



ized) spinors ink{z) and eigenvalues E n k are obtained by 
solving the multi-band Schrodinger equation for the con- 
finement in growth direction of the 2D heterostructure. 
We can then express the spin susceptibility as 



Xij(R,z;R',z') 



d 2 g 
(2^) 2 



Xij(q;z,z') (3a) 



in terms of the 2D Fourier-transformed susceptibility 



.1 



Elk+q - E nk - lfV(] 

where n-p denotes the Fermi function, and 
#^(k,q;z,z') = [^(z)f ■ [j^Ik+qW 

X 6k+q(/)] • [jj^nk(z')\ .(3C) 

Axial symmetry of the 2D system implies E n k = E n k 
and permits the ansatz [23. 24 



ink(z) 



-iJz4>k 



€nk(z) 



(4) 



In 



simplifying calculation of the matrix elements % 
the following, we consider the growth-direction-averaged 
spin susceptibility Xy(q) = J dz J dz' Xy(q;z,z') cal- 
culated at zero temperature. 



Luttinger-model description of quasi-2D holes - Using 
the 4x4 Luttinger Hamiltonian in axial approximation 
for the bulk valence band, the bound states of 2D holes 
confined by a potential V(z) are given in terms of spinor 
wave functions ink(z) that satisfy the Schrodinger equa- 
tion [Jtf + J4?i + Jt? 2 ] ink = E nk ink, where 



,2 r 
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V(z), (5a) 
(5b) 
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dz 



2mo 



71 1 + 71 



K - 4 1 ] - 73 



Ji + Jt 



(5c) 



Here mo is the electron mass in vacuum, hole energies are 
counted as negative from the bulk valence-band edge, and 
we used the abbreviations J± = (J x ±iJ y )/V2, {A, B} = 
(AB + BA)/2. The constants 71 and 7, are materials- 
related band-structure parameters [25[ and depend on 
the quantum-well growth direction. 

Straightforward application of the subband k-dot-p 
method [21, 22 1 yields the energy dispersions E n f. with 
associated eigenspinors ink- At k = 0, the eigenspinors 
are also eigenstates of J z with eigenvalues ±3/2 (heavy 
holes) or ±1/2 (light holes), which are split in energy. 
This phenomenon is often referred to as HH-LH split- 
ting [3, 0|. At finite k, the spinors ink are not eigen- 
states of J z anymore. This phenomenon of HH-LH mix- 
ing arises because a spin-3/2 degree of freedom has a 
much richer structure than the more familiar spin-1/2 
case [2^] . Figure [T] shows the 2D-subband dispersions 
obtained for a symmetric hard- wall confinement of width 
d, using band-structure parameters for GaAs confined in 
[001] direction. The numerically obtained E n k and ink 
serve as input for the calculation of the 2D-hole spin sus- 
ceptibility according to Eqs. ([3]) with Eq. (|4}. 

Results for the 2D-hole spin susceptibility - Using po- 
lar coordinates (q, q ) for the 2D wave vector q and in- 
troducing the scale \o = 2too/(7i^ 2 )j we find that the 
tensor elements of X,j(q) have the generic form 

X xx (q) = Xo [F\\ (q) + G(q) cos (2^ q )] , (6a) 

X yy (q) = Xo [*|| (?) ~ G(q) cos (20 q )] , (6b) 

Xxy(q) =XoG{q) sin(2<?f> q ) , (6c) 

X z M=XoF ± (q) . (6d) 

The functions F^j_(x) and G(x) depend on materials and 
morphological parameters of the 2D hole system, espe- 
cially the hole sheet density n 2 D = ftp/ (27r) , but G(0) = 
generally. Figure [2] shows typical results obtained at low, 
intermediate and high densities where still only states in 
the lowest quasi-2D subband are occupied. 

The density dependence of the 2D-hole spin suscepti- 
bility follows a generic trend. In the limit of low den- 
sity [Fig.^a)], X 2Z (q) > X xx {<\) ~ X yy (<l), and the line 
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FIG. 2. q-dependent spin susceptibility for </> q = and values 
of the 2D-hole Fermi wave vector fcp as indicated in the indi- 
vidual panels. Calculations are based on the band structure 
shown in Fig. [T] 



shape of X zz (ci) is similar to the universal 2D-clcctron 
(Lindhard) result [l3j]. The strong easy-axis response is 
expected [a, H, [t| as a result of HH-LH splitting, which 
favors a spin-3/2 quantization axis perpendicular to the 
2D plane [7| . Interestingly, the behavior of the spin sus- 
ceptibility changes as density is increased. For interme- 
diate values of hole density, results like that shown in 
Fig. EJb) are obtained, exhibiting easy-plane anisotropy 
in the long-wave-length limit and a nontrivial structure 
developing at wave vectors comparable to k-p- The sig- 
nificant deviation from both the univcrsal-2D-electron 



behavior and also the easy-axis response expected from 
HH-LH splitting arises because, as the 2D-hole density 
is increased, higher-A: states get occupied that are more 
strongly influenced by HH-LH mixing. At the highest 
values of density where still only the lowest quasi-2D hole 
subband is populated, easy-axis anisotropy is restored in 
the long-wave-length limit but the response at finite q 
becomes as important in strength as that for q — > 0. 

2D -hole-mediated magnetism - In a magnetically 
doped p-type quantum well, holes can mediate a mag- 
netic interaction between localized moments via the 
RKKY mechanism. The corresponding spin Hamil- 
tonian is given by Eq. ([1]), with Xij{ r a',?p) = 
Xij (R a , z a \ R^, zp) being the spin susceptibility of the 
2D hole system calculated here. For a random but on 
average homogeneous distribution of magnetic ions, stan- 
dard mean-field theory [l2j yields the Curie temperatures 



x,,(q) 



Xo 



(7a) 



q=0 



for ferromagnetic order with magnetization direction par- 
allel to the j axis. The temperature scale 



T n = 



HI 



l) G 2 m 



m 



12 /cb d 7rfi, 2 7i 



(7b) 



depends on the impurity-spin magnitude / and the aver- 
age 3D density ni of magnetic impurities, and its func- 
tional form is that obtained for 2D charge carriers with 
parabolic subband dispersion 0, [2^ corresponding to an 
effective mass too/7i- 

Figure [3] shows the mean- field Curie temperatures for 
perpendicular-to-plane (_L) and in-plane (||) magnetiza- 
tion directions calculated with the same input param- 
eters used for obtaining the subbands given in Fig. [TJ 




FIG. 3. Mean-field Curie temperatures for 2D-hole-mediated 
easy- axis (Tq) and easy-plane (Tj) magnetism as a function 
of the 2D-hole Fermi wave vector &f, obtained for a hard- 
wall confinement with width d and band-structure parameters 
applicable to GaAs. 
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The density dependence of x xx ,zz{^) is directly reflected 

in that of . At the mean- field level, the ordered 
state associated with the maximum transition temper- 
ature will be established. Our results suggest that the 
type of magnetic ordering can be modified by changing 
the 2D-hole density, e.g., by adjusting the gate voltage 
in accumulation-layer devices [28j . Easy-axis magnetism 
prevails at low and high densities, whereas an unexpected 
easy-plane magnetic order emerges at intermediate values 
of the density. Also for high densities, a local maximum 
appears in X xx ((l) at q ps 0.6fcF, which almost reaches the 
value of Xzzi.1 = 0)- Sec Fig. [UJc). Even after averag- 
ing over the polar angle, we find that for kp ~ 1.5 n/d 
the in-plane susceptibility can have a maximum at q ^ 
which is as large as X zz (q = 0). Thus it may be possible 
that the high-density easy-axis ferromagnetic state must 
coexist (or compete) with helical magnetism [12 1. 

Finite-temperature effects - Thermal excitation of spin 
waves (magnons) suppresses the magnitude of the mag- 
netization below its mean- field value Mq. This effect is 
captured by the relation [29|, |30| 



M(T) 
M 



1 



1 



J d\ n q (T) 



(8) 



where nq(T) is the occupation- number distribution func- 
tion of magnon modes at temperature T. A spin-wave- 
rclatcd critical temperature is defined by the condition 



M(T C 



(SW) 



because, for T > T, 



(SW) 



c 



too many 



magnon excitations will have been excited to sustain 
a finite magnetization. In equilibrium, n q (T) is given 
by the Bose-Einstein distribution function ns(£q) = 
l/(e Eq /[ fcBT l — 1), which depends on the spin-wave en- 
ergy dispersion e q . The latter's expression in terms of 
the charge carriers' q-dependent spin susceptibility de- 
pends on the type of magnetic order (Heisenberg, Ising, 
or helical) [29|, but the parameterisation 



(9) 



generally holds for the relevant energy range. The di- 
mensionless quantities eo, v, c v depend on the functional 
form of the q-dependent spin susceptibility and can es- 
sentially be read off Fig. [2J Stability of the magnetic 
order requires both coefficients e~o and c v to be posi- 
tive. Inserting Eq. ([9j) into Eq. ([8]) we obtain an implicit 
equation for the critical temperature. Specializing to our 
situation of 2D-hole-mediated magnetism, we see from 
Fig. dja) that the easy-axis magnetism expected at low 
hole-sheet densities is actually destabilized by magnons 
because £ q < 0. Considering the easy-plane magnetism 
at intermediate densities, Fig. [2jb) reveals that the as- 
sociated magnon dispersion is characterized by v = 2 
and £q = 0, which again implies destabilization of this 
magnetic order due to spin- wave excitations. For the 
easy-axis (Ising) magnet expected at high densities [cf. 




FIG. 4. Real-space spin-susceptibility ^-(R.) of a 2D hole 
system, plotted in units of xo = 27r 2 mo/(7iA 2 d 2 ) as a function 
of kpd for fixed </>r = and kpR — 10. 



Fig. |5fc)], we nn d v = 2 and Eq > 0. In this case a fi- 
nite spin-wave-related critical temperature is obtained. 
Our discussion of finite-temperature effects was based 
on the spin susceptibility obtained for a non-interacting 
2D hole system in an inversion-symmetric quantum well. 
However, recent studies (30| of magnetism mediated by 
2D conduction-band carriers showed that the interplay 
of structural inversion asymmetry and Coulomb interac- 
tions can dramatically increase the stability of ferromag- 
netic order. As the same mechanisms exist for 2D holes, 
more detailed investigation of such effects will be needed 
to conclusively establish the stable types of confined-hole- 
mediated magnetic order. 

Real-space spin response - Fourier transforming Eq. ([5]) 
yields the spin susceptibility in real space Xij i which 
is relevant, e.g., for the proposed use of RKKY interac- 
tions to couple quantum-dot-confined spins as part of 
quantum-information-processing protocols [Til . [Til 17 1. 
Figure [4] shows the density dependence of ^ -(R) for 
a 2D hole system. In the low-density limit, we find 
Xzz ^ Xxx ~ Xyy ~ 0, reflecting the strong easy- axis 
anisotropy arising due to HH-LH splitting As the 

density of 2D holes in the lowest subband increases, in- 
plane components are observed to become relevant. For 
fcp ~ 1.67r/d we find x X x t° be the dominant component 

0. 



uu 



of the spin susceptibility tensor while Xzz ~ X 
Switching of the effective Hamiltonian that describes the 
coupling of two localized spins at a fixed distance is there- 
fore possible by simply adjusting the density of the 2D 
hole system in which they are embedded. 

Conclusions - Our work illustrates the crucial impor- 
tance of valence-band mixing for the spin-related re- 
sponse in confined hole systems, which turns out to 
be much more strongly affected than the density re- 
sponse 31 -33[. It also clarifies the impact of band- 
structure effects on hole-mediated magnetism in 2D sys- 
tems which had previously been only considered for the 
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3D case [3J-(3j| or outside the RKKY limit [37|. 

Acknowledgments - The authors benefited from useful 
discussions with A. H. MacDonald, J. Splettstofier, and 
R. Winkler. 



H. Ohno, |Science 281, 951 (1998)| 

T. Jungwirth, J. Sinova, J. Masek, J. Kucera, and A. H. 
MacDonald, Rev. Mod. Phys. 78, 809 (2006). 
T. Dietl, Nature Mater. 9, 965 (2010). 
K. Pappert, S. Humpfner, C. Gould, J. Wenisch, 
K. Brunner, G. Schmidt, and L. W. Molenkamp, Na- 
ture Phys. 3, 573 (2007). 

A. M. Nazmul, T. Amemiya, Y. Shuto, S. Sugahara, and 

M. Tanaka, Phys. Rev. Lett. 95, 017201 (2005). 

U. Wurstbauer, C. Sliwa, D. Weiss, T. Dietl, and 

W. Wegscheider, Nature Phys. 6, 539 (2010). 

R. Winkler, Spin-Orbit Coupling Effects in Two- 

Dimensional Electron and Hole Systems (Springer, 

Berlin, 2003). 

A. Haury, A. Wasiela, A. Arnoult, J. Cibert, 
S. Tatarenko, T. Dietl, and Y. M. dAubigne, Phys. Rev. 
Lett. 79, 511 (1997). 

T. Dietl, A. Haury, and Y. M. dAubigne, Phys. Rev. B 
55, R3347 (1997). 

S. A. Wolf, D. D. Awschalom, R. A. Buhrmann, 
J. M. Daughton, S. von Molnar, M. L. Roukes, A. Y. 
Chtchelkanova, and D. M. Treger, Science 294, 1488 
(2001). 

D. D. Awschalom and M. E. Flatte, Nature Phys. 3, 153 

(2006) . 

K. Yosida, Theory of Magnetism (Springer, Berlin, 1996). 

G. Giuliani and G. Vignale, Quantum Theory of the Elec- 
tron Liquid (Cambridge U Press, Cambridge, UK, 2005). 

H. Imamura, P. Bruno, and Y. Utsumi, Phys. Rev. B 
69, 121303 (2004). 

M. Pletyukhov and S. Konschuh, Eur. Phys. J. B 60, 29 

(2007) . 



[16] D. F. Mross and H. Johannesson, Phys. Rev. B 80, 
155302 (2009). 

[17] S. Chesi and D. Loss, Phys. Rev. B 82, 165303 (2010). 
[18] J. M. Luttinger, Phys. Rev. 102, 1030 (1956). 
[19] K. Suzuki and J. C. Hensel, Phys. Rev. B 9, 4184 (1974). 
[20] H.-R. Trebin, U. Rossler, and R. Ranvaud, Phys. Rev. 

B 20, 686 (1979). 
[21] D. A. Broido and L. J. Sham, Phys. Rev. B 31, 888 

(1985). 

[22] S.-R. Eric Yang, D. A. Broido, and L. J. Sham, Phys. 

Rev. B 32, 6630 (1985). 
[23] B. Zhu and K. Huang, Phys. Rev. B 36, 8102 (1987). 
[24] B. Zhu, Phys. Rev. B 37, 4689 (1988). 
[25] I. Vurgaftman, J. R. Meyer, and L. R. Ram-Mohan, J. 

Appl. Phys. 89, 5815 (2001). 
[26] R. Winkler, Phys. Rev. B 70, 125301 (2004). 
[27] B. Lee, T. Jungwirth, and A. H. MacDonald, Phys. Rev. 

B 61, 15606 (2000). 
[28] J. C. H. Chen, D. Q. Wang, O. Klochan, A. P. Micol- 

ich, K. D. Gupta, F. Sfigakis, D. A. Ritchie, D. Reuter, 

A. D. Wieck, and A. R. Hamilton, Appl. Phys. Lett. 

100, 052101 (2012). 
[29] P. Simon, B. Braunecker, and D. Loss, Phys. Rev. B 77, 

045108 (2008). 

[30] R. A. Zak, D. L. Maslov, and D. Loss, Phys. Rev. B 85, 
115424 (2012). 

[31] S.-J. Cheng and R. R. Gerhardts, 

|Phys. Rev. B 63, 035314 (2001")] 
[32] V. Lopez-Richard, G. E. Marques, and C. Trallero- 

Giner, J. Appl. Phys. 89, 6400 (2001). 
[33] T. Kernreiter, M. Governale, and U. Ziilicke, 

|New J. Phys 12, 093002 (2010)] 
[34] G. A. Fiete, G. Zarand, B. Janko, P. Redliriski, and C. P. 

Moca, Phys. Rev. B 71, 115202 (2005). 
[35] C. Timm and A. H. MacDonald, Phys. Rev. B 71, 155206 

(2005). 

[36] F. V. Kyrychenko and C. A. Ullrich, 
[37 



|Phys. Rev. B 83, 205206 (2011) 



B. Lee, T. Jungwirth, and A. H. MacDonald, Semicond. 
Sci. Techn. 17, 393 (2002). 



